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A model of the gravitational dipole is proposed in a close analogy to that of the global 
monopole. The physical properties and the range of validity of the model are examined as is the 
motion of test particles in the dipole background. It is found that the metric of the gravitational 
Q\ • dipole describes a curved space-time, so one would expect it to have a more pronounced effect on 

the motion of the test particles than the spinning cosmic string. It is indeed so and in the generic 
case the impact of repulsive centrifugal force results in a motion whose orbits when projected 
on the equatorial plane represent unfolding spirals or hyperbolas. Only in one special case these 
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projections are straight lines, pretty much in a manner observed in the field of the spinning 
cosmic string. Even if open, the orbits are nevertheless bounded in the angular coordinate 9. 
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1 Introduction 



Recently, there has been a flurry of activity in studying exotic cosmic structures, ranging from strings 
to monopoles to walls [1-4]. The motivation to study these objects comes mainly from theories of 
particle physics, where such structures naturally emerge as topological defects which may have arisen 
during phase transitions in the early Universe. Since they are macroscopic, to obtain their complete 
physical description, including global properties, one is inevitably lead to incorporate into this picture 
also gravity^ On the other hand, one can envisage these or similar structures within the general 
relativity (GR) framework only, especially when their field theoretical justification is not available 
or hard to conceive without more sophisticated models, usually invoking additional, hypothetical 
processes. The case at hand is the cosmic spinning string || which can be viewed as a generalization 
of a static cosmic string [1, 2]. Unlike the latter, the spinning cosmic string carries angular momentum 
whose source cannot be convincingly and decidedly established. The only attempt to elucidate the 
origin of the angular momentum made by Mazur || relies on a purely microscopic phenomenon of 
superconductivity. This, however, is not the only problem that one encounters while studying the 
objects under discussion. Some of them are described by sources which themselves are challenged 
by certain pathologies. A good example to call for is that of the global monopole Q. The energy- 
momentum tensor of that object is provided by a respectable theory. However, since the mass of 
the monopole per radial length is constant, its density grows infinitely large when we approach the 
monopole's center. Due to this and the physical requirement that the monopole's mass be finite one 
should think of a realistic gravitational monopole as a shell rather than a ball. 

The still hypothetical objects mentioned in the previous paragraph belong to some sort of extremes: 
the spinning cosmic string in the way it is described by extends to infinity, as does the gravitational 
monopole according to the model of . However, these circumstances should not discourage us from 
considering those models as representing some physical reality to which they are just the first and 
sometimes crude approximation. In fact, it is usually very difficult to find metrics that would give us 
finer descriptions of this reality without the simplicity to be sacrificed, be it a finite cosmic string or 
a global monopole of a finite size. Therefore such models are not unfrequently the most reasonable 
compromises between the complexity that a more accurate description necessarily entails and the 
elucidation of physical contents which one attempts to achieve without pushing the complexity too 
far. Since it is very often the latter that interests us most, it would be rather unwise to discard the 
models that are capable of giving us a glimpse of the contents at a very little expense of means involved. 
As we will see, this also concerns the model we aim at presenting here, the gravitational dipole, which 
shares features of the global monopole both in its infinite spatial extension and a singularity of its 
internal physical characteristics, in this case the density of angular momentum. We do propose to 
overcome these by suggesting some better physical approximation to the strict mathematical model 
in a manner similar to that for the global monopole. 

The nature of spin inducing sources in GR is still a subject of speculations. For example, the Kerr 
solution @ that we will use to derive the metric of the gravitational dipole, believed to be the metric 
of the exterior of the rotating black hole, has not been successfully matched to an interior solution 
that would be generated by some material model, like a spinning shell. Some progress towards this 
goal has been achieved R [TO], though. 



Despite this rather obscure nature of the spinning mechanisms in GR, or perhaps because of this, 



1 The history of topological defects in gravity begins with a seminal paper of A. Staruszkiewicz in Acta Phys. Polon. 
24, 734 (1963), where the defect analogous to the cosmic string is introduced in the 1+2- dimensional context. 
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one is nevertheless challenged to try to understand how the spinning objects can arise and to what 
physical phenomena they may lead in either classical or quantum domain of physics. The goal of the 
present paper is to address the latter issue. In what follows we propose a model of the 3-dimensional 
gravitational dipole in an attempt to investigate the consequences of spinning motions in a general 
relativistic setting. It is in this setting that the gravitational dipole model seems to emphasize the 
role of spin in a more explicit manner than the Kerr black hole. In the limit of vanishing angular 
momentum, the Kerr metric boils down to the Schwarzschild metric, both describing curved manifolds, 
while in the very same limit the gravitational dipole space-time becomes flat. It is in this context 
that the difference between the two spinning space-times is demonstrated most pronouncedly. The 
content of the gravitational dipole metric is purely rotational, without the angular momentum any 
conceivable gravitational interactions are absent here, unlike in the case of the Kerr metric. One 
can therefore hope that through the study of this model one would better understand the physical 
implications of spin in GR in three dimensions which may be different than those due to an essentially 
two dimensional spinning source of the string ||. One would also like to understand how these sources 
can affect the behavior of spin carrying quantum particles. These investigations will be presented in 
another paper [[LI]]. 

The paper is organized as follows. The next section deals with the derivation of the metric for 
the model under consideration and discusses the conditions under which this derivation is valid. It 
also presents the main properties of the model's space-time and considers their physical relevance. As 
explicitly shown in the appendix, unlike in the case of the spinning string the space of the gravitational 
dipole is not even locally flat. Because of this, the spin induced gravity affects the motion of test 
particles in a manner different from the way the spinning cosmic string imparts on them. This motion 
is more complex than in the field of the string as new features characteristic of the centrifugal force 
are present. This is the main observation of the section that follows the derivation of the metric and 
precedes the conclusions that summarize our findings. 



2 The metric of the gravitational dipole and assumptions 
that underline it 

We will derive here the metric of the gravitational dipole and discuss the circumstances under which 
it can be valid as a consistent description of some sort of physical reality. To this end, we start from 
the Kerr metric in the Boyer-Lindquist coordinates || 

ds 2 = §(dt -asm 2 6d(j)) 2 - ^[(r 2 + a 2 )d(j) - ad<P} 2 - § - £W = 
i(A - a 2 sin 2 6)dt 2 + ^f^[a(r 2 + a 2 ) - aA]dtd(f) - %dr 2 - Y 2 d9 2 + (1) 

^[a 2 Asin 2 #- (r 2 + a 2 ) 2 ]d0 2 , 

where S 2 = r 2 + a 2 cos 2 9, A = r 2 — 2Mkr + a 2 , and a — 4r. Here J denotes the angular momentum 
of the rotating black hole, M its mass, and k — ^. We will assume that J = jr, M = mr, j and m 
being some constants and that terms proportional to and of higher order in p — - can be neglected 
in a certain sensible approximation that we want to work out. This amounts to saying that p <C 1 in 
the approximation intended. We will later substantiate this approach in greater detail. We will also 
assume that km < 1, ^ < 1 and because of this the terms kmp (and obviously also kmp 2 ) will not 
be taken into account in a series expansion of the metric (1) in parameters km, kj, and p that seem 
to be the most natural ones for this sort of expansion. 
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Now, from (1) we obtain up to 0(p 2 ) but with 0(kmp 2 ) already discarded that 



g tt — 1 — 2km, — g, rr = 1 + 2km — p 2 sin 2 9, —gee = r 2 (l + p 2 cos 2 i 

= r 2 (l +p 2 ) sin 2 9, -g t4> = 2kjsm 2 9. 



(2) 



This leads us to the metric: 



ds 2 = (1 - 2km){dt - 2kj sin 2 9d<P) 2 - (1 + 2km){\ - E^J)dr 2 - 
r 2 [(l + p 2 cos 2 9)d9 2 + (1 + p 2 ) sin 2 9d<f) 2 } + 0{k 2 j 2 )d<f) 2 . 



(3) 



Neglecting terms 0(k 2 j 2 ) and 0(p 2 ) gives 

ds 2 = (1 - 2km){dt - 2kj sin 2 #d0) 2 - (1 + 2km)dr 2 - r 2 d9 2 - r 2 sin 9 dcp 2 . (4) 

Rescaling (4) by (1 — 2km)~ 1 and changing the coordinates via dp — (1 + 2km)dr brings us to 

ds 2 = (dt - 2A;j sin 2 #d0) 2 - dp 2 - (1 - 2km)p 2 {d9 2 + sin 2 #d0 2 ), (5) 

where we used (1 — 2km)~ 1 = 1 + 2km to be consistent with our assumption km <C 1. Therefore we 
arrived at the metric 

ds 2 = {dt-A sin 2 9d<pf - dp 2 - Bp 2 {d9 2 + sin 2 9d<j) 2 ), (6) 

where A = 2kj and B = 1 — 2 /cm. Since B ^ 1 represents a purely monopole contribution to the 
metric (6) as seen from the metric of the global monopole 0], 



ds 2 = eft 2 - dp 2 - b 2 p 2 (d9 2 + sin 

and we are interested in the purely dipole part of it, we will put B equal 1 in order to focus on the 
dipole content of our derivation. This is tantamount to saying that 

ds 2 = (dt-A sin 2 9d(j)) 2 - dr 2 - r 2 (d9 2 + sin 2 9d(j) 2 ) (7) 

is the metric of the gravitational dipole under the assumptions specified earlier. 

We will now investigate if these assumptions are self-consistent and what is the physical nature 
of the object that produces the gravitational field described by (7). The source of the field is the 
momentum-energy tensor whose components are proportional to the components of the Einstein tensor 
and the latter are worked out in the appendix. Because of their complexity it is rather hard to propose 
a field-theoretical model, or what amounts to the same, some matter Lagrangian that would lead to 
this energy-momentum tensor. This, however, as argued in the introduction, should not be considered 
a major obstruction. The metric (7), due to the assumption that the angular momentum per unit 
radial length is constant, describes an object whose total angular momentum is infinite. To make 
this situation look more physical one can think of (7) as a good approximation to the metric of an 
exterior of a shell of thickness d in its vicinity, a distance r 3> a from the center, which is required 
to satisfy the condition p< 1. Then the angular momentum is finite and proportional to d. We will 
soon show that under reasonable and physically natural circumstances one can make a = 1 and since 
making a = 1 coincides with imposing on r the condition r > 7-^171 the discussed condition p < 1 
is not restrictive at all. As a matter of fact, because km <C 1, p < 1 is easily met. Why it should 
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be a shell rather than a ball can be justified by equilibrium considerations. For the equilibrium to 
take place the centripetal force present in any rotating system should not overcome the gravity force, 
that is for a test particle of mass /i there should be fxu 2 (r)r < kM ^ ^ , where uj{r) and M(r) are the 
angular velocity and the mass of the matter that generate our metric. This means that 



9 / \ JvJ\ / f 9 / \ 

uj (r) < — — = kmr . (8) 



On the other hand, from the definition of the angular momentum J = J r 2 u(r) dM(r) one obtains 

— — = 7 = const = r u)(r)— — = r u(r)m (9) 
dr v ' dr v ' w 

according to the assumptions employed. Eq. (9) requires io(r) ~ \, thus leading to a sigularity with 
r — > 0. It is the main reason why we want a shell as our model for the metric of the gravitational 
dipole. The situation here is somewhat similar to the mass distribution in the global monopole, where 
being constant implies p(r) ~ which holds also for the model of the gravitational dipole under 
study. Obviously, u cannot grow to infinity without destroying the stability of the whole configuration. 
It is Eq. (8) that demands for uj{r) — S% a minimum radius r such that ^ = ^ or ro = \/lm' ^ 
choosing (7=1, which can always be accomplished by rescaling the length, we see that the balance 
of the configuration is maintained as long as its inner radius is greater than ^ fc J^ 1/2 . This choice of C 

leads via (9) to a = — = 1. 

A few remarks on physical properties of the introduced space-time seems to be in order here. First 
of all, the metric (7) as an axially-symmetric one possesses two Killing vectors d t and 8$. In the 
region r < AsinO, the latter, which has closed orbits, becomes timelike. Any observer following these 
orbits would experience causality violation. Fortunately, in the regime in which the metric of the 
gravitational dipole is valid (that is km <C 1 and so also kj 1 for a = — = 1) this region lies inside 
the shell discussed above and its physical relevance is only apparent. As long as one stays outside of 
this shell, one can dismiss this issue. 

As shown in |~2| the spinning cosmic string indroduces the time delay between the arrival time 



of two particles moving in opposite directions that follow an arbitrary closed trajectory around the 
string. One can rightly expect a similar effect in the case under study. The only difference between 
these two physical situations is that in the metric of the gravitational dipole this delay depends on 
the angle 9. For light rays describing circular paths with = Qt, r = R =const, and 9 =const, the 
eikonal equation ds 2 = implies Q + = fj ,_ A 1 sin ^ e and Q- = — R+J ^ sin 2 e with the Q + assigned to the 
ray moving along with the direction of rotation. Assuming R > Asm 2 9, which, as argued earlier, is 
always valid for our model, one obtains for the time delay 



T = 2ir 



1 1 



10 J ft_ 



(10) 



The last comment concerns the very name we chose for the model presented. As is well-known, the 
gravitational charges are only alike, thus preventing the existence of "electric" gravitational dipoles. 
However, the dipole we introduced in this paper is of "magnetic" nature. Since in GR every form 
of energy contributes to the gravitational field, in the case under consideration it is the rotational 
( "magnetic" ) energy of the dipole that does so. As opposed to the spinning cosmic string that carries 
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a nonzero flux of the "magnetic" field, but does not produce a detectable field itself f], the gravitational 
dipole creates such a field as manifested by the curvature of its space-time. 



3 The motion of test particles 

As observed in the previous section, the metric of the gravitational dipole possesses two Killing vectors, 
d t and d<p that can be used to find integrals of motion of test particles in the space-time of the dipole. 
The integrals themselves are useful in simplifying the equations of motion for the particles by casting 
them into a dynamical system. To employ this strategy let us recall that 

(d a \0 = const = —P a (1) 

along the geodesic £, where d a is a Killing vector, £ is a vector tangent to the geodesic, and P a is a 
momentum conjugate to the coordinate x a , chosen so that d a is tangent to it. Eq. (1) is valid also 
for a zero-mass particles, i.e., when /i = except for its last equality. The brackets in (1) denote the 
scalar product of two vectors separated by a vertical dash. 
Eq. (1) when applied to an axially-symmetric metric yields 

E 

gut + gt<j>4> = n = ( 2a ) 

Hoc 

9tti + 9<jxp<P = — = A, (26) 

fJLQC 

where 7 and A are constants related to the energy E and the azimuthal angular momentum I of the 
test particle. When used for our metric Eqs. (2) lead to 

XA — 7X 

t = ~2 > ( 3fl ) 

= ^, (36) 

where X = r 2 — A 2 sin 2 9. The remaining equations of motion can be obtained from the Hamilton- 
Jacobi equation 

^1^-^ = 0. (4) 

In this formalism P M = J^- and since we have already established that P t = —^ocy and = /iocA 
one easily arrives at 

/ A 2 ^ 2 e \ 1 P 2_ P 2_2A 1 



and then at 



^ + ^L + 7 2 ^ 2 sin 2 ^+r 2 
HqC z sm I 



+ ( 7 2 - 1) + —j 1 



(6) 



2 For this reason the spinning cosmic string could also be called the gravitational fluxon. 
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By introducing a separation constant K such that 



J ^ + -4r+7 2 ^sin^ = K 2 (7) 

one obtains 

# = (7 " 1) + — — • (8) 

With the substitution x a = ^ Eqs. (7) and (8) can be reduced to the following dynamical system 
for r and 9 



/ X 2 \ 1/2 

0= K 2 -^- 7 2 A 2 sin 2 ^ . (96) 
\ sin 9 J 

Making the substitution z = cos9, this last equation can be transformed into the energy con- 
servation equation for the one-dimensional motion of unit mass particle in the potential V(z), or 
explicitly 

£ + V(z) = E (10) 

with E = \[K 2 - A 2 - 7 2 A 2 ) and V{z) = az 2 + bz 4 , where a = \[K 2 - 2^ 2 A 2 ) and b = \~i 2 A 2 . One 
more substitution, x = z 2 , brings us to the equation 

x 2 = 8x(E - ax - bx 2 ) (11) 

which when supplemented with the constraints x 2 > and x > describes the physical motion in the 
coordinate 9. 

Depending on the values of E and a, there are different ranges of 9 in which this motion takes 
place. Let us first consider the case of a > 0. As can be easily seen, the motion is possible only for 
E > 0, E = representing a singular situation with the orbit constrained to the equatorial plane 
9 — | . For E < both X- and x + are negative and upon consulting with (12) one sees that the motion 
is impossible (x 2 < 0). When the "total energy" E is positive, the trajectory is bounded between 
this plane and the plane determined by the equation Q_ = arccos ^fxZ, where x. = For 

2 

a < there are more options and only if E < E min = — |r no motion takes place. For E = E min 
there exists a stationary orbit in the plane 9 = arccos with x 2 = This corresponds to the 

previous situation of E = 0. When the "total energy" exceeds E min there exists an orbit bounded 
by the planes 9_ = arccos y/xZ, and 9 + = arccos \J max[0, x + ] , where x + = — a+%/ ° 2 & + — . This also 
includes the case E = for which 9 + — |. 

If both a and E are negative the solution to (11) can be expressed in terms of elliptic functions. 
To see this let us rewrite this equation as 

x 2 = — 8bx(x — x+)(x — X-), (12) 

where x± are defined above. Upon the introduction of k 2 = — and y such that the last 

equation becomes 

y 2 = -2bx4l-y 2 )(l-k 2 y 2 ), (13) 
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Through the substitution u = i^/2bx + t one can bring (13) to the standard form of the differential 
equation for elliptic functions 

whose solution is y(u) = sn(u + 5) with 5 being some constant. Finally, the solution to (9b) is found 

to be 

9{t) = arccos y/x+ sn ( J2bx-it + 5) ■ (14a) 



If E = and a < then x + = and (12) reduces to 



' dy s 
du , 



2bx + y 2 (l-y 2 ), 



where y 2 



Upon the substitutions u = ^j2bx~Zt and 1 — y = w the last equation boils down to 



solution to which is given by w = ±tanh(u + c), c being some constant. Eventually we find that 



9(t) = arccos 



1 - tanh 2 (c + ^2bx_t^j 



(146) 



For E positive and independently of a, the trajectory is bounded between the equatorial plane and 
the plane 6_ = arccos a -^ a ^+ 4bE ^j _ This case can be reduced to the previous one by changing 
x + , which is always negative, to — \x + \, and replacing x by y via x = y — \x + \. In doing so one ends 
up with 

y 2 = -8by(y-(x_ + \x+\)) 

that looks exactly like the equation leading to (14b), except that x_ gets replaced by x_ + \x+\. 
Obviously, the solution to (9b) for such conditions is 



6(t) = arccos y (x_ + ^1 — tanh 2 + \j2b(x^ + |rr + |)t^ — \x + \ 



(14c) 



However, the fact that the motion is bounded within some range of 9 does not necessarily mean 
that its orbits are closed. As a matter of fact, we will now demonstrate that the orbits of test particles 
are open. To this end let us find the trajectory of the test particle as a function r ((/>). By combining 
(3b) and (9a) one obtains 

d * = r*lp?-Plr*y (15) 

where C = A — 7 A, a 2 = 7 2 — 1 > 0, and f3 2 = K 2 — 2A\ r y. Let us now consider three cases of f3 2 > 0, 
(3 2 = 0, and f3 2 < 0. They will lead to different classes of trajectories. In the first case, changing 
variables \ = u subsequently u — ^ cos x brings one to the solution 



a 



r(0) \(3\ 



cos[\(3\(<f> + B)/C], 



(16) 
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where B is a constant of integration. The case /3 2 = is the easiest to work out. Through only one 
change of variables £ = u, it immediately leads to 

«*> = (17) 

In the last case, the change w = p sinh x enables one to arrive upon integration at 

' =-^smh[-\P\(<i> + B)/C\. (18) 



r(0) 

As seen from (3b), C > ensues that increases with time. However, the sign of C has no bearing 
on the character of motion, all it does is to change its direction. The trajectories in the last two cases 
describe unfolding spirals which asymptotically approach infinity in the limit — > B. For C > this 
is achieved through increasing values of that can run from — oo while for C < the orbit develops 
towards decreasing values of <fi with <fi running from +oo. That it is so is determined by the condition 
f{4>) > 0- 111 the first case the trajectories are hyperbolas unless (3 = ±C which corresponds to E — 
and results in a motion along a straight line in the equatorial plane. Both the hyperbolas and the 
straight line make the closest approach to the center a distance ^ from it. The asymptotes of the 
hyperbolas are defined by the equation |/3|(0 + B)/C = ±|. When the motion is projected onto the 
equatorial plane, the asymptotes divide this plane into four sectors with the particle moving in the 
one that satisfies the condition cos [\(3\((f) + B)/C] > 0. 

The equations of motion for particles of zero rest mass can be worked out in a similar manner. In 
fact, if the mass term is dropped in (5) the only change it brings about is in Eqs. (6), (8), and (9a), 
where 7 2 — 1 is replaced by 7 2 which represents the square of the total energy of the particle. The 
way this change affects the presented solutions is really insignificant. 



4 Conclusions 

We presented the model of the gravitational dipole in a yet another attempt to investigate the impact 
of spinning sources on gravity in a general relativistic framework. The model studied has some 
advantages in this respect over the Kerr black hole as its gravitational field is of pure rotational origin, 
thus allowing one to examine this problem in its disentanglement from the main gravity source, the 
mass. We showed that the model has sensible properties if it is thought of as an approximation to a 
rotating shell. The space-time of the dipole is curved, therefore it should not come as a surprise that 
the way gravity associated with the curvature affects the motion of test particles is more pronounced 
and complex than that of the spinning cosmic string. Unlike for the spinning defect-free string, the 
trajectories of test particles in the field of the dipole when projected on the equatorial plane are not 
straight lines anymore, the only exception being the motion in the plane itself. All trajectories are 
open and bounded in the 6 coordinate. 
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Appendix 



The only non-zero covariant components of the Einstein tensor for the metric of the gravitational 
dipole are as follows 

3a 2 cos 2 9 
G tt - — 4 , 

a 2 cos 2 9 



Gaa — 



a 2 cos 2 9 



a 2 (5r 2 cos 2 9 - 3r 2 cos 4 9 + 3a 2 cos 2 9 - Qa 2 cos 4 9 + 3a 2 cos 6 9) 
G 4>4> ~ ~a ■ 

a(r 2 — r 2 cos 2 —3a 2 cos 2 9 + 3a 2 cos 4 9) 

G<t>t = Gt<t> = — A • 

The Ricci scalar R is 

2a 2 cos 2 9 
K = . 

r 4 

Clearly, the space-time of the gravitational dipole is curved. 
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